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■ Abstract 

A weak-coupling scaling diagram for the Lyapunov exponent and the integrated density 
I of states near a band edge of a random Jacobi matrix is obtained. The analysis is based 

■ on the use of a Fokker-Planck operator describing the drift-diffusion of the Priifer phases. 

(N ■ 
> 

o 

O : 1 Main result and short overview 
o 

-p^ . This work considers one-dimensional discrete random Schrodinger operators, typically given by 
p ^! a periodic background operator and a weakly-coupled random potential. The periodic (non- 
random) operator has a band structure. In the vicinity of a band edge, we present a rigorous 
pertubation theory in the coupling constant and the energy for the Lyapunov exponent and 
the density of states. This leads to a new and rich scaling diagram describing these two self- 
averaging quantities at a band edge. 
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The scaling behavior at energies within a band has been well understood for a long time. 



' For sake of concreteness, let us focus on the Anderson model H = A + W given as the sum of 
^ the discrete Laplacian A and a disordered potential V = Xlnez"^" coupled with a (small) 

coupling constant A > 0. Thouless [Thoj found a perturbative formula for the Lyapunov 
exponent (inverse localization length) at an energy E = 2 cos(/c) in the band spectrum of A: 

8 sm (k) 

where Va- is one of the identically distributed and centered entries of V and E^r denotes the 
expectation value. A rigorous proof was provided by Pastur and Figotin using the Priifer 
phases, also called the Dyson-Schmidt variables |PF] . This amounts to justifying the so-called 
random phase approximation, stating that the Priifer phases are distributed according to the 
Lebesgue measure on the unit circle. A further generalization was given in [JSSj . One feature 
of ([1]) are the singularities at the band edges A; = 0, vr of the unperturbed operator A. In fact, 
the control on the error terms in ([1]) breaks down at those energies, as does the random phase 
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approximation. The latter also happens at the band center, namely for k = f, this leads 
to anomalies in the perturbative formula ([1]) first found by Kappus and Wegner [KW] . and 
consecutively analyzed by several authors [DGl IBKl ICKt ISVW] . A perturbative formula as ([1]) 
for the band center with a control on the error terms was proven only recently |Slj . 

This paper considers the band edges and shows how the is modified, leading to the 
new scaling diagram near a band edge. Parts of this diagram were already given by Derrida 
and Gardner [DGJ . These authors actually found the correct scaling in the parabolic regime 
of Theorem [1] below, but the wrong prefactor (c/. the comment at the end of Section [8]). 
Moreover, they could not give a better justification of their scaling Ansatz than that it leads to 
a differential equation they could solve. Our more conceptual approach shows why the scaling 
is natural in the situation considered in |DGj . It exhibits a far richer scaling behavior near a 
band edge and also allows to rigorously control the higher order corrections. Moreover, if one 
wants to give a perturbative proof of uniform positivity of the Lyapunov exponent in an energy 
interval around a band edge, all of the scaling behaviors considered below are needed. 

Let us now describe the main result in more detail. We consider a one-parameter family 
of random Jacobi matrices {H\^i^)x>q given as the sum of an L-periodic background operator 
Ho = Hq^w and a random perturbation H\^^^ — Hq which is linear in the (small) coupling constant 
A. More precisely, for every fixed configuration uo, the operator Hxj^ acts on G £^(Z) as 

Hx,u.\n) = tx,u.{n + l)\n+l) +Vx,M\n) +tx,u,{n)\n-l) , nGZ, (2) 

where the coefficients tx^^^{n) > and v\^^{n) G M are constructed as described in the following: 
let (ti, . . . , ti, Vi, . . . ,vl) be given real constants with > for / = 1, . . . , L; set S = [—1, 1]^^ so 
that each a G S is of the form a = . . . , iL{a),Vi{a), . . . , VL^a)); then f2 = x {1, . . . , L} 

is the configuration space and to each u = ((crm)mez, k) E Q there are associated sequences 

tx,w{k - 1 + Lm + I) = ii + Xii{am), vx,uj{k - I + Lm + I) = vi + Xvi{(7m) , 

which for A sufficiently small satisfy t\^^{n) > 0; these sequences define H\,^ by ([2]). In order 
to make {Hx,uj)x>o into a family of random operators, we equip Q with a probability measure 
P = X Y2t=i where p is a probability measure on S. Expectation values w.r.t. P and p 
will be denoted by E and We suppose that Eo-(ti(cr)) = Eifj{vi{a)) = 0. 

For every fixed energy E E W and coupling parameter A, there are two self-averaging 
quantities of interest, namely the integrated density of states (IDS) Mx{E) and the Lyapunov 
exponent (or inverse localization length) 'yx{E). The definitions will be recalled in Section H] 
below. The periodic operator Hq = Hq,^ has a band structure and we are interested in the 
scahng of the IDS and Lyapunov exponent at one of its band edges Ef, (band touching excluded). 
In order to state the precise result, we need to introduce two quantities. Let be the random 
transfer matrix at energy E over a unit cell of length L (see Section [2] for the explicit formula). 
Then set 

X = dETTi%%)\^^^^ , = 9.Tr(r5)U=o. 
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Theorem 1 Let Eb be a band edge {band touching excluded) of the periodic background operator 
Ho = Hq^w of a family of random Jacobi matrices {H\^i^)x>q. The perturbation is supposed to 
be non-trivial in the sense that x^ does not vanish p-almost surely. Then the scaling near the 
band edge is 

Xx{Eb + eX^) = ^^oiEh) + + 0(A°+^) , ^xiE^ + eX^^) = B + 0(A^+'^) , (3) 

where A, B, a, (3 and 5 > depend on e G M and rj > as described in the following and resumed 
in Figure 1. 

(i) (Elliptic regime) Let ?7 < | and let the sign of e ^ be such that E^ + eX^ is inside the band 
of Hq. For the case of the Lyapunov exponent, we also suppose r]> ^. Then 

V 

a = -, (3 = 2-r] , 
ii) (Parabolic regime) Let rj = ^, then 



and 



ex\ 1 
A = sgn(e) , B = — 



E^(|x^ 






\ex\ 





a = [3 
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and A = A{e) and B = B{t) are given by integrals written out explicitly in Section [HI For 
rj > ^, A and B are independent of t, namely the result is the same as for e = 0. 

(iii) (Hyperbolic regime) Let | < 77 < | and let the sign of e be such that Efy + e X^ is 
outside the band of Hq. Then 



and 



a > - , p = - 
2 ' ^ 2 



B = ^/\^\ 



The result is summarized in Figure [H All the error estimates, in particular, the value of 5, 
are controlled more explicitely and given in Section [HI As will become more clear below, the 
terms elliptic, parabolic and hyperbolic regimes reflect the nature of the corresponding averaged 
transfer matrix to lowest order in A. Let us point out that we do not provide the asymptotics 
for the IDS in the hyperbolic regime, but only give an upper bound on it. Actually the IDS 
between the Lifshitz tails and the band edge of the periodic operator is exponentially small in 
A, a fact that results from a more delicate large deviation behavior. This will be dealt with 
elsewhere. The Lifshitz tails themselves were analyzed in |S2j . 

The proof of Theorem [1] begins (in Section [2]) with an adequate basis change on the transfer 
matrices. A band edge is characterized by the fact that the modulus of the trace of the transfer 
matrix is 2. If there is no band touching, the transfer matrix is then not diagonalizable and its 
normal form is hence a Jordan block (Appendix A). The action of a Jordan block on projective 
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space (in fact, identified with the Priifer phases) has exactly one fixed point, which is unstable 
under random pertubations induced by the randomness in the Hamiltonian. In order to analyse 
the associated invariant distribution on projective space, a further basis change blowing up 
the vicinity of the fixed point is necessary. If this is done adequately, a formal perturbative 
calculation typically leads to a differential equation for the invariant distribution. 

Instead of studying the concrete form of the random matrices obtained in Section [2] after the 
various basis changes, we rather choose to single out the more general concept of an anomaly 
of a family of random matrices ( cf. Section [3]) . The term anomaly is chosen by reference to 
the center of band anomaly in the one-dimensional Anderson model as studied by Kappus and 
Wegner |KW] . The latter is a special case which is analyzed in a prior work |S1] . The formalism 
to calculate the Lyapunov exponents and rotation numbers at a general anomaly is built up in 
Section |H In particular, it turns out that one needs to evaluate certain Birkhoff sums for the 
perturbative calculation of these quantities. The perturbative evaluation of these Birkhoff sums 
directly leads to a differential operator on projective space (Propositions [1] and [2]). In some 
situations this operator is of first order, and consequently the anomaly is called first order as 
well (this is dealt with in Section[5]). However, in the more interesting cases treated in Sections[6] 
and [71 the differential operator is of second order and of the Fokker-Planck type. The latter 
operator was already used in [ST] in order to study the center of band anomaly, but here one 
is confronted with the supplementary difficulty that its ellipticity is destroyed at a band edge. 
Thus a thorough analysis of the corresponding singularities is needed. Section |6] deals with the 
Fokker-Planck operator and its groundstate, while Section [7] shows how it is used to calculate the 
Birkhoff sums. Both sections heavily depend on Appendix B, where inhomogeneous singular 
first order ordinary differential equations are studied in detail. Sections [3] up to [7] are kept 
slightly more general than needed for the proof of Theorem [T] as completed in Section [HI but 
we hope that this stresses the structural aspects of the analysis. There may possibly also be 
further applications. 

Acknowledgment: We are thankful to Andreas Knauf for numerous comments and discus- 
sions. This work was support by the DFG. 
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2 Normal form of transfer matrix at a band edge 



In this section we motivate and carry out the basis changes on the transfer matrices in the 
vicinity of a band edge. The final normal forms obtained in the various regimes will then 
motivate the definition of the term anomaly in the next section. 

First, one rewrites the Schrodinger equation Hx^^ip = Eip, with G M, G £^(Z) and 
i^A.o; as given in ([2]), in the standard way using transfer matrices {e.g. |JSS] ). Due to the 
periodicity of Hq it is convenient to introduce the transfer matrix over L sites associated to one 
configuration cr G S of the disorder on these sites: 

= n ^""(^'^ + >^i,{a)A + , r^M = (^^ ) . (4) 

1=1 ^ ^ 

Note that T^{t,v) G SL(2,R) and hence T^^ G SL(2,M). If Eb is a band edge of Hq, then 
Tr(7^^^) = ±2 and, because there is no band touching , OeTtIT^^^) ^ 0. One can show that 
necessarily Tff^ has only one eigenvector (see Appendix A). Hence one can find a basis change 
N G SL(2, M) such that NT^^^N"^ is a Jordan block with eigenvalue either 1 or —1. This 
elementary fact was already efficiently used for the analysis of the Lifshitz tails |S2j . Using this 
basis change, one can write the full energy rescaled transfer matrix as follows: 

N r^r'' = ± exp + J2 A" Qv.. ) • (5) 

Here the exponents rjk are in N + r^N, are put into increasing order rj^ < rjk+i and the Q,,^.,^ 
are in the Lie algebra sl(2,M). Hence the lowest terms are X'^Qr],a or XQi^^, pending on the 
value of T]. Let us note that Qkri,cT are independent of a unless kr] G N. Also, Ei^{Qi „.) = 
unless there is a such that krj = 1. Each combination of the signs can occur, the one inside 
the exponential indicates if it is a lower (+ sign) or an upper (- sign) bandedge. For sake of 
concreteness, we choose both signs positive. 

A few further remarks on the Qr]k,cr will be relevant later on. Because there is no band 
touching at E^, for e 7^ one gets 

^ edETTiN%lN~')\^^^^ = Tr(^Q,,. (^J ^ = ex , (6) 

where the first identity follows from the definition ([5]), Duhamel's formula and the cyclicity of 
the trace, and the second one defines x. By the same calculation (with e = 0) 

dxTTiNT,';^N-')\x=o = Tr(^gi,. (^J |^ ^ = , 

the latter by definition of the centered random variable x^-. Of course this expression is centered 
because the perturbation is centered. However, the assumption in Theorem [1] on the non- 
invariance of the band edge under perturbation translates into 

2^ 



E 



Tr Qi,^ 



1 1 
1 



> . (7) 



It turns out (see Section H] or [JSS] ) that one needs to study the random dynamical system 
given by random iteration of the natural action of the transfer matrices NT^^^^^"^ on 
S]^ = M/vrZ, given by (see Section H] for details) 



± 



X.(TrL Pn — 1 

I A.fJn t^n — 1 I 



cos(e) 
sin(^) 



At A = 0, the action is given by the graph of Figure El It has an unstable parabolic fixed 



Figure 2: Dynamics at A = at lower band edge, i.e. + sign in 



point at = 0, which is being approached from the left under the dynamics without passing it 
though. For A > 0, the dynamics is the same to lowest order A°, but the random perturbations 
of higher order may allow the dynamics to cross through 6* = 0. Nevertheless, the angles 9^ 
are close to ^ = for most n and thus their distribution is more and more concentrated in the 
vicinity of 6* = as A tends to and, in fact, converges weakly to a Dirac peak. In order to 
extract the shape of the non-trivial distribution it is necessary to rescale the neighborhood of 
6' = with an adequate power of A. This is done by conjugation with 



A,<5 



A'' 
1 



The conjugation has the following effect on 2 x 2 matrices: 



N: 



\,5 



a b 
c d 



A 



-5 





c 



+ 



a 
d 



+ X' 



b 




(9) 



One deduces from ([5]) 



exp 



N: 



A,<5 



^ ^ k>l 
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Note that the powers of A appearing in the exponential now he in N + ?7N + {—6,0,6}, but 
none is negative for 6 < min{l,r]}. The foUowing three contributions of low order in A turn 
out to be relevant: 

Xs,i = \' J) , Xs,2 = X'^-' ]imX'Ny,iQ,,.N-,\, Xs,.,3 = \'-' hmX' Ny,,Q,,, N-,\ , 

where the notation reflects that Qj,,^ is independent of a. By equations ([9]) and IQ, Xs^2 has 
a non- vanishing constant entry only in the lower left corner unless e = 0. Similarly by and 
d?]), has a centered random entry in the lower left corner with positive variance, while 

all the other entries vanish. It is helpful to think of Xs^i and Xs^2 as drift terms in the action 
(IH]), and of X5^o-,3 as generator of a diffusion. Now we blow up the parabolic fixed point by 
augmenting 6 until the drift X5 1 is of same order of magnitude as either the drift Xs^2 (then 
6 = 2) or the diffusive force (variance) of Xs^a,3 (then 2(1 — 5) = 6 so that 5 = |). Hence we 
choose 

= mm <—,->. 
12' 3/ 

Explicitely, for 77 < |, one has with this choice of 6 

N.,NT-:;-^^'^N-^N^} = exp (aI (^°^ j) + A^i + Q'^ + AQ',. + , 

(10) 

where O = 0{X~, A^+2) and the specific form of Q'^ and Q[^^ will not be of importance for 
the asymptotics of Theorem [TJ For rj > 1, the second term is of lower order than the first, but 
as it is centered, we shall see that it does not enter into the asymptotics of the IDS and the 
Lyapunov exponent because its variance is smaller than (^(Aa). Hence (11 01) is an anomaly of 
first order in the sense of the following section. In the case ^7 = 1, 

iV„iVT5->-'iVll ^ exp(Ai [I i) . At [^^ J) . 0(A)) . (11) 

This is actually the same formula as (llOp . but it was written out again in order to clearly 
exhibit the so-called anomaly of second order (see Definition [2] below), namely the variance of 
the first centered term is precisely of the same order of magnitude as the expectation value of 
the second. If 77 > |, one obtains the same formula with e = 0, but with an additional error 
term 0(A''"i). 



3 Definition of anomalies 

Let us consider families (Tx^a) \m,(Tes of matrices in SL(2,]R) depending on a random variable 
a in some probability space (S, p) as well as a real coupling parameter A. In order to avoid 
technicalities, we suppose that Tx^a has compact support for small A. Furthermore we assume 
that the dependence on A can be expanded in some power series with non-negative exponents 
(not necessarily integers). Later on, we shall choose Tx^a to be the matrices in (fTOj) and (ITT]) . 
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Definition 1 The value X = is an anomaly of the family (Ta,o-)agm,(tge if for all a E and 
a sign that may depend on a E T,: 

To,. = ± 1 . (12) 

In order to further classify the anomahes, let us suppose that for a fixed A- and cr-independent 
basis change M e SL(2,R) to be chosen later, the transfer matrix is of the following form: 

MTx,M'' = ± exp ( 5^A'"^ P^„. J . (13) 

\A.>1 / 

Here Pr,^,a ^ sl(2,]R), r^^ > for A; G N and K, such that rjj < rjk if j < k, rjx = 2rii, 
TjK+i < + and p(P^j.,o- = 0) < 1 for A; = 1, . . . , i^T — 1 (which means none of these matrices 
shall be identically for p-almost all a). 

Definition 2 Let (TA,o-)AeM,o-es have an anomaly and suppose given the expansion ( |T3l) for 
a fixed basis change M. The anomaly is said to be of first order and ICth kind if K, < K , 
E(P^^,cr) = for k = 1, . . . ,/C — 1 and E(P^^,o-) is non-vanishing. An anomaly of first order 
andKth kind is called elliptic i/det(E(P^^,cr)) > 0, hyperbolic z/det(E(P^^,o-)) < and parabolic 
z/ det(E(P^^,o-)) = 0. Note that all these notions are independent of the choice of M . 

If E(P^j,,cr) = 0, for k = 1, . . . , K — 1 {then the variance of P^^-^^^ is non-vanishing) , then an 
anomaly is said to be of second order. 

It may happen that the transfer matrices have to be modified as follows in order to obtain 
an anomaly. If there is a mapping A 1— > M\ G GL(2,]R) for A > where M\ is independent of 
a, such that 

limMATA,.M,-i = ±1, (14) 

then we say that Tx^a is transformed to an anomaly by Ma. Note that the limits limAj^o ^a and 
limAjo need not exist. 

As argued in the last section, the transfer matrices at a band edge can be transformed 
into anomalies flTU]) and flTT]) of respectively first and second order. The term anomaly first 
appeared in the work of Kappus and Wegner |KWj on the band center of a one-dimensional 
Anderson model. Indeed, as discussed in |Slj . the square of the transfer matrix at a band 
edge is an anomaly in the above sense (of second order, if the random potential is centered). 
More generally, higher powers of transfer matrices are transformed into anomalies for rational 
quasimomenta, and thus lead to anomalies in higher order perturbation theory [DG]. A more 
systematic classification of anomalies as in Definition [2] was done in [Slj . However, we felt 
it adequate to change the term degree of an anomaly used in [Si] to the present order of an 
anomaly because it leads to a differential equation of corresponding order. 

4 Phase shift dynamics 

The bijective action St of a matrix T G SL(2, R) on Si = M/vrZ = [0, vr) is given by 
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± 



Teg 
WTeeW 



cos{e) 



(15) 



with an adequate choice of the sign. This defines a group action, namely Stt' = StSt'- In 
particular the map St is invertible and S^^ = St-'^- Note that this is an action on projective 
space. In order to shorten notations, we write 

At an anomaly one has 5a,o-(^) = + 

Given an initial angle 6q and iterating this dynamics defines a Markov process ^^(cu, ^g) on 
= X {1, . . . , L}, i.e. for oj = ((cr„)„g2, k) one defines iteratively 



(16) 



In order to shorten notations, we denote 6'„(a;,6'o) by 9n- Under adequate identifications, this 
corresponds to the dynamics ([H]). The 6n are also called modified Priifer variables. 

In order to analyze the dynamics in more detail, some further notations are needed. For 
A; e N, we define the trigonometric polynomials: 



Pk,a{0) 



53m 





p 






ee) J 



1 

7! 



1 



More explicitly, in terms of the matrix elements one obtains 
n.,. = ^T'^) =^ Pk,M = c,„.cos2(^)-^„.sin2(e)-a,„.sin(2^). (17) 

Now starting from the identity 



\MTx,M-'\ee) 



{v\MTx,M-'\ee) ' 
the definition f[T^ and the identity {v\eg) = e*^, one can verify that 



\U {v\ee) 



As one readily verifies that 

Pl^, = -det(P,„,)l, 53m 
it follows that 



2{v\ee) 



{v 






ee? 


2{v 





+ C(A^^'+i) . 



(t;|ee)2 



-Pi,adePi,a{0) 



^ 1 

Sx,M = + $^A*pfc,.(e) + -X^'^'pi^adePiAO) + 0(A^^+O • 



k=l 



9 



As one has exp A***^ -Pr^^.o-) ^ = exp (— Yl P'nk,^) same procedure leads to 

Sxl{(^) = ^ - Y.^''VkAO) + i^\'''Vnu<^deVr„,M + 0{\^-^^) . (19) 

k=l 

Our main interest is the perturbative calculation of the Lyapunov exponent 7(A) and ro- 
tation number 7^(A) associated to the random family of matrices {T\^an)n>i which are defined 
by 

N-l 

7(A) = lim -E^E.log(||MT,,.M-Vll), (20) 

n=0 

and 

1 1 ^ 

7^(A) = - hm -Ej2^,^,A0n) . (21) 

vr N—*oo iV 

n=l 

where ipx,a : ^ R is the phase shift given by ^x,a{0) = SxA^) - ^- If ^a,^ = NT^;^'^" N-^ 
where the latter is the transfer matrix as given in (jl]) associated to a random Jacobi matrix 
H\^uj with an L-periodic background operator defined in ([2]) and N is an arbitrary basis change 
such as the ones used in ( ITOl) and (fTTI) . then the Lyapunov exponent 7a (-E) (inverse localization 
length) and IDS M\{E) of the random Jacobi matrix are 

7A(i5 + eA") = ^ 7(A) , Ux{E + eA") = - i 7^(A) mod ^ . (22) 

The first identity follows immediately from the definition; the second is a consequence of the 
oscillation theorem and the gap labelling for the periodic operator (see e.g. [JSSj for details). 
In both of the expressions fl20l) and fl2T|) . one can now expand each summand w.r.t. A. 

Lemma 1 Set 

as well as ^ 

Then p^,j{9) = 53m(a^^o- — /3^o-e^*^). Furthermore, with errors of order 0{X^^^), 

k>l 

(23) 

+ %e ('\'?fci+'7fc2 (3 '3 + A e^'^ -3 8 e^'^)^ 

and up to errors of order (9(A''^"'"^^) 

= ^ra 5^(a^'= - P,,,.e'^') + i A^'"^ (-2 a,,,. e^^^ + Pl^e''^')) . (24) 

fc>i 
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Proof. This follows from straightforward algebra using 

{ee\T\ee) = ^Tr(T) + 3?e ( {v\T\v) e^'') , Tr(|P,,.|2 + P^j = 4|^^^^|2 ^ 

as well as the identities Tr(P^^o-) = and (■u|P^,^|t') = 0. □ 

Once these formulas are replaced in ( 120|) and ( l2Ti) . one hence needs to consider Birkhoff 
sums of the type 

/^(/) = T7E5^/(ft.), /(/) = lim J^(/), (25) 

n=0 

for TT-periodic functions /. For 7(A) and 7^(A) one actually only needs the functions e^*^ and 
g4ie_ ^}^ggg Birkhoff sums have to be evaluated perturbatively in A with a rigorous control on 
the error terms. For this purpose, one needs to know the distribution of the 6n as generated by 
the random dynamics (fTSjl . This is the main focus of the next sections. 



5 First order anomalies 

It turns out that it is easiest to calculate the Birkhoff sums lN{f) in case of an elliptic first 
order anomaly of /Cth kind. As then det(E(P^^^o-)) > 0, one can choose the basis change M in 
( IT3i) such that 

E(P....) = (° "0^) , (26) 

with /i 7^ 0. Hence E(p/c,cr(^)) = Z^- 

Proposition 1 Let Tx^a have an elliptic anomaly of first order and ICth kind and suppose that 
E(pk:,(7(^)) = A*- Then for any f e C^{Sl), one has 

Mf) = r^/(^^) + 0(A'"c'-^^(A'^-iV)-l) , 
Jo ^ 

with IC' = min{}C, K} > K, where K, is such that P^, is the first uncentered term in (fT3!) after 



p 



Proof. Because InU) = c+Iiq{f —c) for c = JJ" dO f{6)/Tx, we may assume that JJ" AO f{6) = 0. 
Then / has an antiderivative F G C^{Sl). Using the Taylor expansion and 2?7i = rjx, 

( K \ K-X 

V fc=l / k=l 

From this one deduces after taking expectation and summing over n: 

In{F) = Jjv(P) + /iA'"^ lAf) + 0(A'"=',iV-i) . 
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As /i 7^ this proves the estimate. 



□ 



Concerning hyperbohc and parabohc first order anomahes, we do not provide an exhaustive 
treatment, but rather present some procedures on how these anomahes can (possibly) be trans- 
formed into a more accessible anomaly such as one of second order. This apphes, in particular, 
to the hyperbolic first order anomaly (fTOj) with ex negative, corresponding to the hyperbolic 
regime in Theorem [H For a hyperbolic first order anomaly of /Cth kind, let us first choose M 
such that 

'-/i 0' 
/i. 



(27) 



Then to lowest order in A, = is the unstable and 6 = ^ the stable fixed point of (the 
averaged dynamics of) iSa,<j- Like in (fTTj) . we use the A-dependent basis change N\^s- Let x be 
the smallest exponent such that E(c^^^) > 0. First suppose that E(cg^o-) = for all ^ < ^ri/c + x, 
then choose S = x — hvic > so that 2(x — 5) = rjic- One obtains 



A,. 



± exp 



A2 



VIC 







+ ... + A 



VK 



-a 



+ 



{2i 



where, moreover, E(6, 



and E(a 



VK,(rJ 



—fi < 0. This is a second order anomaly for 
which the Birkhoff sums are analyzed in Theorem H] below. If, on the other hand, there is a 
^ < ^Vk. + X such that E(c^ g.) ^ and ^ is the smallest such exponent, choose 6 = — rj/c > 
so that 



N- 



A,<5 



± exp 







... + A 



VK 



'VK,<^ "VK~S,<^ 



-a 



'VK,<^ 



where ( is either r/i or x — 6. For both cases we have 2^ > rjic- As E(^^_5o-) = 0, one thus 
has again a hyperbolic first degree anomaly. Repeating the procedure, one may now be in the 
above advantageous case. 

A parabolic first order anomaly of /Cth kind can be treated just as a the band edge (which 
is parabolic to 0th order) in Section [2J One first chooses M to be a rotation matrix, such that 
the parabolic fixed point is 6* = 0, namely E(P^^^o-)eo = 0, and then considers the matrices 
Nx^s MTx,aM~^ N^j. Let again x be the smallest exponent such that E(c^_^) > and let ^ 

for all 



be the smallest exponent such that E(c^_o-) 7^ 0. If ^ > |x + Ivic and E(a 
1] < + Vic), then one chooses 6 = \{2x — Vic) > in order to get a second order anomaly 



ri,a ) 



± exp 



y^iix+VK) 



x-<5,o- 



c 







+ 



6 Fokker-Planck operator of a second order anomaly 

In order to calculate /iv(/) at a second order anomaly, we first refine the proof of Proposition [H 
This naturally leads to an associated differential operator, as shown in the next proposition. 
The remainder of this section analyzes the properties of these operators. 
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Proposition 2 Let the random family {T\^cr)cr€T. have a second order anomaly and F G C^{Sl). 
Then, for 

f = Eipl^)F" + (2E(p^,.) + E(pi,.9,pi,.))F' , (29) 

one has 

Ij^(f) = C>(A™,(A''^Ar)-i) _ 
with fj = min{?7j^, r/i + 772} where K is the first uncentered term in f|T3|) after Pq^^a- 
Proof. A Taylor expansion implies 
/ k-i 

y k=i 



fc=i 



Taking expectation values and summing over n, one gets the estimate by the same argument 
as in Proposition [1] because rjx = ^rji. □ 

Let us sketch the further strategy. In view of ( l29l) . it is natural to introduce: 



p = E{pI^) , q = 2EipK,a) + Eipi^^depi,^) . (30) 

Both p and q are the trigonometric polynomials of degree 4 on S*^, namely linear combinations 
of e^^*^, e^^*^ and the constant function. Moreover, p > 0. Hence p can have at most one zero 
of order 4, or two zeros of order 2. Actually a zero 9 of order 4 only appears in the rather 
special case described next. By applying a rotation matrix M in ( fT3l) . one may assume 9 = ^ 
and then one readily checks that 

piO) = E(4,J cos^(^^) , P,„. ^ ° ° 



This is precisely the situation in ffTTj) . Furthermore, let us note that q can have zeros with 
orders adding up to at most 4. 

Proposition [2] now states that one can control BirkhofF sums of functions which are in the 
range of the differential operator on C^{Sl) given by 

C = {pde + q)de . 

The formal adjoint on L^(S'^), also considered as an operator defined on C^(S'^), is given by 

C* = de {de p - q) ■ 

The operator C* is a Fokker-Planck (or forward Kolmogorov) operator on S\ describing a drift- 
diff'usion dynamics on S"^, while C is the associated backward Kolmogorov operator [Ris] . The 



13 



operator C* was already used in |Slj . but only in the situation of a strictly positive p. In the 
case of the anomaly (ITTl) corresponding to the parabolic regime at a band edge, p has a zero. 
Hence the ellipticity of C* is destroyed. This is the main technical problem that has to be dealt 
with in the present work. 

As Ran(£) cKer(£*)-'-, a necessary condition for controlling the Birkhoff sum lN{f) is 
hence that / GKer (£*)-*-. As will shortly be shown in Theorem [21 in many situations Ker(/}*) 
is one-dimensional and spanned by a smooth positive, L^-normalized function p. Let us call 
p the groundstate of the Fokker-Planck operator, even though the spectrum of the Fokker- 
Planck operator is non-positive. Furthermore, the necessary condition / G Ker(£*)-'- combined 
with / G C^{Sl) will in most cases turn out (Theorem [3] below) to be sufficient for finding a 
solution F G C^{Sl) of the inhomogeneous differential equation fl29|) . Thus one deduces from 
Proposition [2] that, for / G C^iSl), 

iNif) = J d9pi9)fi9) + 0(A^-+-^-,(A''-iV)-^) . 

In some other situations, relevant for the hyperbolic regime in Theorem [1], the kernel of C* is 
given by a Dirac peak p (if C* is considered as an operator on the space of distributions), and 
the same formula holds. 

Before proving the main result on the groundstate p, let us point out that p can be seen 
as the lowest order formal approximation to the Furstenberg measure on 5"^ associated to 
family (Tx^a)a£j: of random matrices. It is defined by p3Lj 

r ux{d9) f{9) = E r ux{d9) fiSxAO)) > / ^ C{Sl) . (31) 
Jo Jo 

In fact, supposing h'xid9) = px{9) d9, (131!) leads to 

E{deS-^l^{9)px{S~^],m) = PxiO) ■ 

Furthermore assuming that px is twice differentiable and px = p + 0{X'^"^+^~'^'^), this gives the 
equation C*p = (c/. [Si] for further details). 

Theorem 2 The Fokker-Planck operator C* has a unique groundstate p, which is non-negative, 
normalized and continuous except possibly in the following cases: 

(I) p has a zero 9 of order 2 {and possibly another zero) for which, moreover, 

dep{9) = q{9) = 0, d^gp{9) > deq{9). (32) 

(II) p has a zero 9 of order 4 for which, moreover, 

E{pkA0)) = , E{de PK,.m < . (33) 

Furthermore, the groundstate is infinitely differentiable except possibly at one point 9, at which 
p has a zero of order 2 and 

dep{9) = q{9) , d^,p{9) < deq{9) . (34) 
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Let us note that (1321) is equivalent to 

Furthermore, let us point out that whenever there is a point 6 which is a zero of both p and 
q, then the Dirac peak 6q is a groundstate of C* if considered as an operator on distributions. 
This Dirac peak is dynamically stable under the drift-diffusion described by C* if, moreover, 

d'iidep) - q){e) > 0, I e {1,3}. 

This is precisely the case in (I) and (II). A corresponding dual result, namely that the Birkhoff 
sums /iv(/) are equal to f{6) to lowest order, will be proven in Theorem H] below. Let us 
also note that there may be two such angles {e.g. if E{pK,cr) = 0), in which case the stable 
groundstate is degenerate. Finally, we encourage the reader to draw the curves of p and q in 
each of the cases of the proof below, and interprete p as a diffusive force, and g as a drift. This 
allows to better understand the formal proof, as well as the result itself. 

Proof of Theorem [2l The equation C*p = can be integrated once. Hence one needs to solve 

[pdg + q) p = C , q = (dgp) - q, (35) 

where the real constant C has to be chosen such that (1351) admits a non-negative, vr-periodic, L^- 
normalized solution p. The equation fl35l) is locally precisely the type of singular ODE studied 
in Appendix B, with coefficients which are moreover real analytic. The only supplementary 
property of fl35l) is that the inhomogeneity is constant. When writing out solutions, we use 
similar notations as in the appendix, except that we shall put tildes on w and W in order 
to distinguish them from the corresponding dual objects appearing in the proof of Theorem [3] 
below, and that we do not include the constant inhomogeneity C in the definition of W below 
as done in (j47j). 

If p > so that there is no singularity and C* is elliptic, the groundstate p can readily be 
calculated as in [Si]. Set, for some 9 E Si, 

^m^j(d^^, (36) 

Then 

p = C,e-^ {C2W + I) , C^= ^, C = CiC2, (37) 

W[p + vr) 

where one first determines the normalization constant C\ which then also fixes C in (!35|) by 
the last identity. Note that W{6) > for > so that, in particular, is well-defined. 
Furthermore C2 = if w is vr-periodic. The groundstate p is clearly real analytic in this case. 

Now let d be the only zero of p. It is either of order 2 or 4. If g(^) 7^ 0, one is locally near Q 
in the situation of Proposition [3t^iv) or (v). For each C in (1351) . the smooth integrable solution 
is unique to one side of d. Continuing it cyclically around S\ shows that it is also unique on 
the other side of Q. For 9 < 9 < 9 + it, this smooth solution is given by 

p = Ce-^W , (38) 
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with w and W defined as in (l36ll using ^<6'<^ + 7rin the first equation and 9 = doi6 = d + Tx 
in the second one, such that the singularities cancel. It only remains to choose C 7^ (and 
hence the ODE (!35l) itself) such that p is normalized. 

If p has two zeros §1 and 62 of order 2 for which q{6i) 7^ and 9(^2) 7^ 0, one has to 
distinguish two cases. If the signs of q{6i) and 5(^2) are the same, then one proceeds cyclically 
twice as in the case with one zero of p which was just treated. The solution is as in fl38l) and 
C 7^ 0. If on the other hand Q < 61 < 62 < n with q{6i) < and 5(6^2 ) > 0, then one chooses 
C = and sets 

= k<0<k, ,39) 

'^^ ' \ otherwise , ^ ' 

with w defined as in (136|) for some 6 with §1 < 9 < §2 and some normalization constant C3. 
Note that this p has no singularities because hm^j^^^ w{9) = +00 and lim^^^^ w{9) = + oc. As 

shown in Proposition [3](iv) and (v), p is smooth at both §1 and ^2- 

Now follow the singular cases where p and q have common zeros. The situations (I) and (II) 
are of this type, but there will appear others for which we then show as stated that there is a 
unique continuous groundstate. As argued in Appendix B, there can then only be integrable 
solutions of ( l35l) if the equation is homogeneous, namely C = 0. One such solution is p = 0, but 
there may be others. First let p have only one zero 9 of order 2 and suppose that q{9) = 0. Then 
^ is a zero of order 2 of g if and only if dBq{9) = 0, that is E(9g p\^^(9)) = 4 E{de PK,a{9)). In this 
situation, there is a one-parameter family of smooth solutions according to Proposition [3](i). 
However, these solutions are typically not periodic and hence this is part of situation (I). If 9 
is a zero of order 1 of g, then Proposition [3](ii) and (iii) can be used for analyzing the local 
regularity. Case (iii) applies if dgq{9) < 0. Thus one has a two-parameter family of continuous 
solutions near 9; however, one parameter has to assure the periodicity of the solution so that 
it remains continuous, while the second has to be chosen so that p is normalized. Hence 
the unique continuous groundstate is p = C^e'^ with an adequate normalization constant 
C4. In the situation (I) when ( l32ll holds with a strict inequality, one has dgq{9) > so that 
Proposition [3](ii) implies that the solution is unique. As p = is one solution, it is the only 
one. 

Next let us consider the case where p has two zeros 9i and 6*2 of order 2, and that at least 
one of them, say 9i, is a zero of q. If deqipi) > (namely situation (I) again), there is only the 
zero solution due to Proposition [3](ii). 

If doqipi) = and ^(^2) 7^ 0, the singularity 9i can due to Proposition [31^i) only be resolved 
by choosing C = 0, but then the zero ^2 for which q{92) 7^ leads to a singularity. Thus there 
is no continuous solution and this is included in situation (I). If deq{9i) < and say q{92) > 0, 
then one can appeal to Proposition [3](iii) at 9i and Proposition [3](iv) at 6*2, and construct the 
groundstate exactly as in fl5I?]) . However, p might not be differentiable at 9i. If 9i and 6*2 are 
zeros of both p and q, then either deq{9i) > or dgq{92) > 0. Indeed, q = ^ dgp — 2E{pK,a) 
and Ei{pK,a) is a trigonometric polynomial of degree 1 and hence can only compensate either 
dep{9i) > or dep{92) > 0. Therefore this corresponds again to situation (I) and there is no 
normalizable solution. 
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In the last remaining case ^ is a zero of p of order 4 (hence the only zero of p) and q{6) = 0. 
If deqiO) < 0, which is equivalent to E(c?6i PK,a{0)) > 0, Proposition [3](vii) can be applied and 
the groundstate p = is smooth. If dgq{9) > 0, then Proposition Myi) implies that the 

zero solution is the only solution. If deq{0) = and E(px,o-) = identically, then q = j dgp has 
a zero at 9 of order 3 and dgq{9) > 0, so that Proposition [3](ii) implies that the zero solution is 
the only solution. The same holds if dgq{9) = and Ei{dgPK,cT{9)) ^ due to Proposition [3](iv) 
or (v), apphed either to the left or right of 9 pending on the sign of E((9gPA',o-(^))- The last 
three cases are regrouped in situation (II). □ 



7 Birkhoff sums at second order anomalies 

In this section we calculate the Birkhoff sums lN{f) perturbatively for many second order 
anomalies. This involves the analysis of the operator C dual to the Fokker-Planck operator. 
We will not present a treatment as exhaustive as Theorem [2] for C* and do not treat the singular 
cases where p and q have common zeros, except for the case of Theorem H] below, which is needed 
for the analysis of a band edge. Comments on the remaining cases are given at the end of the 
section. 

Theorem 3 Let {Tx^^)^^^^, have an anomaly of second order. Suppose q{9) ^ whenever 
E(Pio-(^)) = 0- Furthermore let p be the groundstate of the Fokker-Planck operator given by 
Theorem [21 Then, for f e C^{Sl) one has 



TT 



Mf) = / d0/(^^)p(^^) + 0(A™,(A'^-iV)-i) 







with 7] = mm{rif^,rii + 772} where K is the first uncentered term in f|T3l) after Pqj^^a- U^iPia) 
has no zero, it is enough to suppose f G C^{Sl). 

Proof. As already argued in Section [6], one can control the Birkhoff sums using Proposition [2] 
only if / GKer(£*)-'-. Under the hypothesis stated, the kernel of C* is one-dimensional by 
Theorem [2] and spanned by the groundstate p. Replacing / by / — (p|/), we may now assume 
that Jq d9 f{9) p{9) = and have to show for such / that the Birkhoff sum Iwif) is of the order 
stated. This will follow directly from Proposition [2] once we have found G G C^(S'^) satisfying 

{pde + q)G = f , r d9G{9) = 0. (40) 

7o 



Indeed, the second identity allows to take an antiderivative F of G which then satisfies ( 129|) . 

First let us consider the case where p has no zero and C 7^ in (l35il . Setting, for some 
9 G Si 
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the solution is given by 



G = e-"" {W + ci) , ci 



-w(e+Tr) 



where e~"'(^+'^) 7^ 1 as this happens if and only if w and w are vr-periodic leading to C = 0. One 
readily checks that the first equation in (BOl) is satisfied and that G is periodic, hence in C^{Sl) 
because / G G^{Sl). Furthermore, multiplying the first equation in fHOj) by p and integrating 
over Si shows that also the second equation holds: 

= J pf = J p{pde + q)G = - J G {dep - q)p = -C j G{e) . 
Next, if p has no zeros and C = in fl35|) . then C2 = in fl37j) and 



p = Cie-'^ = Ci — , {pde + g) e-'" = , (42) 
p 

so that W defined as in (jH]) satisfies W{d + n) = W{6) =Q because / pf = 0. Hence 

G = e-"- {W + C2) 

is a (periodic) function in (7^(5^) and, for adequate choice of C2, of vanishing integral. 

Now let p have one zero 6. The hypothesis implies q{6) 7^ 0. According to Proposition [3](iv) 



or (v), the first equation in (HOl) has a C^-solution in a neighborhood of 6, which is unique to 
one side of 9. Continuing this solution around 5*^ one deduces that the solution is unique, and 
actually given by G = e~^W. Due to the argument in the proof of Theorem [2l C 7^ in fl35l) 
in this case, so that by the same argument as above the second equation in fj^Oj) is satisfied. 
The case with two zeros 61 and 62 of p for which q{6i) and 5(^2) have the same sign is treated 
similarly ( cf. the proof of Theorem [2]) . 

Finally we deal with the case of two zeros < ^1 < ^2 < of p for which q{Oi) > and 
9(^2) < 0. By ( l39|) . the support of the groundstate is [6*1,^^2] and hence not all of Si- For 
9 G [6*1, 6*2], we choose 9i < 9 < 62 in the first equation of fHT!) and 6* = 6*1 in the second 
equation of (HT]) . Then G = e~^W is a solution in (^1,^2)- Due to Propostion [3](iv) this is 
the unique solution that can be smoothly extended to the left of 9i. Furthermore (142!) and 
f pf = imply W{92) = 0, so that this solution can also be smoothly extended through 6*2 
by Propostion [2]^v). On both sides one has one free parameter. One is chosen such that G 
is vr-periodic and the remaining one such that the second equation of fHOl) is also satisfied. 
This reflects that Go = ^~^X[e2 §1] ^ smooth positive solution of the homogeneous equation 
{pde + q)Go = 0. □ 

The remaining singular cases for which Theorem [2] guarantees existence of a unique ground- 
state cannot be treated by the techniques of Theorem [3] without further hypothesis on /, which 
seem somewhat unnatural, but reflect the delicate dynamical behavior at such points {e.g. one 
needs that / and its derivative have a particular behavior near the common zero of p and q). 
The next theorem deals with an anomaly corresponding to situation (II) of Theorem [21 albeit 
with a strict inequality in (133|) . 
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Theorem 4 Let {Tx^a-)a£T: have an anomaly of second order. Suppose that 9 is a zero of order 
4 of E(p^^) and a zero of order 1 of 'Ei{pKa) o,nd that, moreover, E{dePK a{9)) < 0. Then for 
f e C'isi) 

iNif) = fie) + 0(A™,(A^-iV)-^) 
with fj = min{?7^, r^i + 172} where K is the first uncentered term in ( fT3l) after Pr)K,cr- 

Proof. The procedure is exactly as in the proof of Theorem [3l Hence we search for a solution 
of (140!) for a function / satisfying f{9) = 0. Close to 9 the equation {pdg + q)G = f has a two- 
parameter family of solutions due to Proposition [3]^ vii). One parameter is fixed by requiring G 
to be periodic. The other refiects that one can, by the same procedure, find a non-vanishing 
solution Go to [pde + q)GQ = 0. As Go has definite sign, one can use it to normalize G so that 
also the second equation in (l40l) is satisfied. □ 

As to extensions of Theorem HJ the case with an equality in the second equation of (15^ 
cannot be treated by the presented technique because {pdo + q)Go = has only the trivial 
solution by Proposition [3](iv) or (v). The case with one zero of order 2 corresponding to 
situation (I) in Theorem[2]can be treated similarly as long as 3 E((9g pf^ai^)) < 4 |E((?6i PK,a{d))\, 
a condition which originates in Proposition [3]^iii). 



8 Application to a band edge 



This section contains the proof of Theorem [H Hence let Ta,o- = Nx^rN T^^'^'^^^ N ^N^ l be the 
anomaly given in (fTOj) or (fTTj) . Then (!22l) allows to calculate the Lyapunov exponent and the 
IDS, so that one may focus on the calculation of 7(A) and 7?.(A) based on (120|) . (12T!) and the 
expansions given in Lemma [TJ 

Let us begin with the elliptic first order regime, hence (fTOj) . The adequate basis change 
assuring (12^ is 

'1 

,0 — |e x| 



M 



Then {MNx,sN)T^^^'^ {MNx,sN)-^ is equal to 

exp (x^exl'^ (0 "j) - A^"^-^ (° °) + A^P, + AP,. + 0(A^,A^+i; 



Here E(Pi_H,a.) = 0, and also E(Pi o-) = 0. With the notations of Lemma [H this gives Pr = 

and /5i-|,o- = «Xo-|ea;|~5/2. Thus Proposition [1] implies /(e^-'*^) = (9(A2,A^~^) for j = 1,2. 
Hence the first 4 terms with powers A^* in (1251) either vanish or are of higher order, the same 
holds for the mixed terms A^'=i"'"^*2, the first of the terms with A^'"= is also of higher order, but 
the second one gives the leading order contribution: 



7(A) = A^-^^ + 0(A^,A^-^) 
8 ex 
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which together with (122|) estabhshes the first claim. Similarly one verifies = «|ex|2 so that 
by ([21D one finds 7^(A) = A? |ex|^ + 0(A^ A^"''). 

In the hyperbolic regime, the basis change leading to (1271) is 



M = 





— \ex[ 


(1 





Then {MNxrN)T^^^^'^\MNxrN)-^ calculated from ([TO]) is equal to 

exp (^A^|ex|i ("^ + A^"i^ ("^ "/) + ^"^^ + + (^(A-,A^+i) j . 

Now let us apply a further basis change A^a,<5 as described in fl28l) . namely in the present situation 
Vic = 2^ X = ^ ~ 2 = 1 ~ Due to ([9]) the transformed anomaly then becomes of second 

order: 

exp(Af|4„|4 0) , Ai|..|l ("^ «) ^ A-f , A-? . o) , 

where O = 0{X^ , A^~^). Here Pi^r „ and -Pg-Sa possibly of lower order than the second 

term (but not as the first one as r/ < |), but they are both centered so that they do not enter 
into the asymptotics (c/. Definition [2]) . The obtained second anomaly can now be dealt with 
by Theorem m implying that /(/) = /(|) + C(A"?~\A^^^). Furthermore, one calculates 
(3ii^tj = —iXa\ex\~^ /A, 72,0- = xl\ex\^^/8 and Pr = —\ex\^. Replacing into fl23|) one realizes 
that the factor of A^4 vanishes so that 

7(A) = A2 |ex|5 + 0{x'-^-\X^--i) . 

Concerning the rotation number, one has an = and I 2/4. The latter implies 

that the contribution A^* vanishes in (12^ due to the ^m. Thus 7^(A) = C(At'-\ Ai-5). Note 
that the condition ^ — 1 > | is equivalent to > |. 

Finally let us deal with the parabolic case for which ffTTl) is directly a second order anomaly. 
The polynomials (130|) are explicitly given by 

p{e) = E{xl) cos\e) , q{e) = ea; - l + (ex + l) cos(2e) -2E{xl) cos^O) sin(0) . 

Because g(|) = — 2 7^ 0, Theorem[2]guarantees the existence of a unique normalized groundstate 
p G C°°{Sl). As ai^^ = iXcr/2, (3i^^ = —iXa/2 and 71 = xl/2, one deduces from Lemma [H 
and Theorem [3] that up to errors of order 0(A) 

7(A) = xl f^{ex + 1) I dep{e) sin(2 0) + j dO p{e) {I + 2cos(2 0) + cos(4^))') , 



20 



and, furthermore that 7?.(A) is, also up to errors of order 0{\), equal to 
^ (^ex - 1 + {ex + 1) j de p{e) cos(2 0) - j dO p{e) {2 sm{2 9) + sm(40))^ . 

Let us remark that this formula is not the same as equation (35) of [DGj . even though the 
invariant measure of [DGj coincides with the above. In fact, their expansion of the expression 
log(l + A^/^t) ~ A^/^t is erroneous for large t. 



Appendix A: Transfer matrix at a band edge 

Let us consider a real analytic family E ^'Rh^ G SL(2, M) with 

Tr(r°) = 2, 9£;Tr(r^)U=o 7^ 0. 

If is the transfer matrix of a periodic Jacobi matrix, this corresponds to a band edge at 
E = Q. We show that necessarily the eigenvalue 1 of T° has geometric multiplicity 1, so that 
the Jordan form of T° is non-trivial. A similar statement holds if Tr(T°) = —2. For this 
purpose, let us choose the notation 

^ - [c + CE d + DEj + ^(^)- 

Then one deduces ad — be = 1, a + d = 2 and A + D 7^ 0. Furthermore the order E of 
det(T^) = 1 implies that aD + dA- cB -hC = We need to show that 7^ 1. If a 7^ 1, this 
is true. Hence suppose a = 1 so that also d = 1. Then ad — he = 1 implies that either 6 = or 
c = 0. Hence aD + dA - eB -hC = Q implies that either A + D - cB = ox A + D -hC = 
Because A + D ^ Q it follows that either c 7^ or 6 7^ 0. 



Appendix B: Inhomogeneous singular first order ODE 

This appendix recollects some results about the solutions of the inhomogeneous first order 
ordinary differential equation (ODE) for a function y = y{x) on the interval U = (a, b) 

py' + qy = r , (43) 

where p,q,r G C"^{U) and m G N U {oo}. In particular, we are interested in smooth solutions 
of the singular case where p, q and r have zeros at x G ?7 of order < lp,lq < m and 
respectively. We suppose that p and q have no further zeros in U. If Ir < min{/p, Ig}, then there 
cannot be any solution y G C^{U) because then 

-y + -y = ^ 

leads to the contradiction = 1 in the limit x ^ x. Hence let us suppose that l,- > min{/p, Iq}. 
In the regular case Ip < Iq, the homogeneous equation py' + qy = has a one parameter family 
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of solutions e where w{x) = | G C™ ^^'^^{U) is any antiderivative. Then a one-parameter 
family of solutions of fj43|) is obtained by the method of variation of constants: 

y(x) = e-"'(") r ds e"'^^) 44 , (44) 

where the integral is well-defined because Ir > Ip- In the singular case, one has Iq < Ip and 
Iq < Ir and the zero of p effectively decouples the intervals (a, x) and {x, b). Any antiderivative 
w{x) = j p diverges as x x. Nevertheless, we shall show below that these divergences may 
cancel out in (jH]). Let us note right away though that any solution y G C^{U) has to satisfy 

y(x) = lim ^ . (45) 

This is a boundary condition for both intervals (a, x) and (x, 6). The question is whether there 
are differentiable solutions (at the point f , in particular) and whether such a solution is unique, 
or there are one-parameter or two-parameter families of such solutions. 

Proposition 3 Let Ir > mm{lq,lp} and set I = Ip — Iq. Then the ODE ( H3ll has continuous 
solutions y with a number of free parameters as described in the following: 

(i) -(f^ < 0, there is a one-parameter family of solutions y G C"^~^^^^{U). 

(ii) If I = 1 and 9|(x) > 0, the solution y G C"^~'"'{U) is unique. 

(iii) If I = 1 and d^{x) < 0, then there is a two-parameter family of solutions y G C"{U) as 
long as n < m — Iq and n < \d^{x)\^^ . 

(iv) If I > 2 is even and d^^{x) > 0, the solution y G C™^'«(t/) is unique for x > x, but there 
is one free parameter for x < x. 

(v) If I > 2 is even and <9'|(x) < 0, the solution y G C™^'''(t/) is unique for x < x, but there 
is one free parameter for x > x. 

(vi) If I > 3 is odd and 9'|(x) > 0, the solution y G C"^^'"'{U) is unique. 



vii) If I > 3 is odd and 9'|(a;) < 0, there is a two-parameter family of solutions y G '"'{U). 
Any other solution y has a non-integrable singularity at x 



Proof. In case (i) the solutions are given in (jH]). Next let us consider the case (iv). Dividing 
( H3!) by q gives ^y' + y = ^. Hence it is sufficient to consider the case py' + y = r with 
p, r G C''^^'"'{U) and Ip = I. Let w be any antiderivative of ^ on f/ \ {x}. Because of d^p{x) > 
and / is even, one has \im.x]xw{x) = -\-oo and lim^j£w(x) = — oo. Moreover, the growth of 
1 1/7(0;) I is at least I because I > 2. Thus for all n G N 

lim p"(x)e"'(^) = + 00 , lim = . (46) 

xTx xlx p[x) 

Now e""^ is a solution of the homogeneous equation py' + ?/ = 0. Because \mixixw{x) = —00, 
the only possible continuous solution oi py' + y = r for x > x is 

y[x) = e-"'^^) W{x) , W{x) = [ ds e""^'^ 44 ' (47) 



p{s) ' 



22 



where the last integral exists due to (H6|) . For x < x, let W be any antiderivative of e""^ and 
also set y = e~^W . Hence y is now a solution of py' + y = r for x ^ x. It remains to be shown 
that y G C'^~'"'{U). First of all, y is continuous because 

W(x) e"'^^) ^ 

lini = lini = lini . . ^^"^ = r(x) = , 

x^x x^x ' x^x '—r-T 

p{x) 

the latter by (H5|) . In the second equality, de FHospital's rule could be applied because for x | x 
the numerator and denominator both converge to by (jl2D and fH7|l . while for x | x as the 
denominator converges to oo by fHBl) . 

Next follows an inductive argument in order to check the continuity of the higher derivatives 
of y at X. Let us set: 

Qn = np' + 1 , r„ = r^_i - q'^_iy^'''^^ , Tq = r . 

Note that g„ has no zero at x. One can check by induction 

+ gnZ/(") = r,. (48) 

Let ty„ = w+n log(|p|). Now e"""" is a solution of the homogeneous equation py^"'^^^ +qn y^"-* = 
and satisfies due to fj46l) 

lim p(x) e"'"^^) = + oo , lim e"'"^^) = . (49) 

x'lx xlx 

Now let y*^"'^ n < m — Ig, he continuous by induction hypothesis. Then r„ G C^(?7). Set 
Wn = y^"^ e"*". Due to (jl9]), one has limaj^^ W^„(x) = 0. Using the identities 

(n+l) ^n^ " ~ QnWn , , ^ 

which follow from (H5|) and = — , one obtains 



g«(a^) W^n(a^) + rn{x) e"'"(^) 
" ■ ■ x— >£ e"'"(^'*p(x) 

g^(x)Vr„(x) + r^(x)e"'"(^) -^U^) l/^''^(^) + ^n(^) 



lim ?/^"^^''(x) = lim 
lim 



«„(x)(p/(a;) + g„(x)) gn+i(^) 



where de I'Hospital's rule can be used for x J, x, because the numerator and denominator 
converge both to according to fj49l) and for x t 3; as the denominator converges to ±oo. 
Therefore y^"^^ is continuously differentiable in x. 

The case (v) is identical, once the sides x > x and x < x are exchanged in the above 
argument. In case (vi), one has lima.^£ ti;(x) = — oo (namely, for both x J, x and x t x), so 
that one has to construct the solution as in fH7|) both for x > x and x < x; hence the solution 
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is unique. In case (vii), one has lima;^^ = + oo and the homogeneous solutions on both 
sides vanish at x faster than any power, so that one has a two-parameter family of solutions. 

Finally let us consider / = 1 and first the case (ii). Hence lim^.^^ i/7(x) = — oo so that 
lim^.^^ gW'Ca;) — g_ Now the second limit in fH6|) may diverge to + oo, but one readily checks that 
^ is integrable so that W can be defined as in (H7!) for x > x and x < x. Now the continuity 
is proven as above, and higher derivatives can also be treated as above because 

lim e'""^^) = 0, 

SO that both nominator and numerator of ?/("+^) converge to and de I'Hopital's rule may be 
applied as before. 

In case (iii), lim^.^£w(a;) = oo so that lima,^^. e~"'*^^^ = 0. Furthermore ^ is never integrable 
and the two-parameter family of solutions y is defined as in dH]). The main point to note is 
that Wn = {p'{x)~^ + n) log(|x — x\) + 0{1), hence 

lim e"'"-'^''^ = oo , lim p(a;) e""""'^^^ = - oo , 

x'\x xlx 

under the hypothesis stated in (iii). Thus de I'Hopital's rule can be invoked to calculate 
lima;^£ 1/*^") as above (with the index n shifted by 1). Let us remark that in case |p'(a;)|^"'^ G N, 
one can prove better differentiability, but only for a one-parameter family of solutions (just as 
for the ODE xy' = Ny, N e N). 

The singularity of all other solutions is at least of the type \x — in cases (ii) to (vii). □ 
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